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THREE-DIMENSIONAL  PLASTIC  PLOW  UNDER  UNIFORM  STRESS* 

William  Prager 

Brown  Universityf  Providence,  R.  I#,  U,  S«A* 


Abstract*  In  view  of  the  role  that  regions  of  uniform 
stress  play  in  tne  theory  of  plane  plastic  flow,  the  most  general 
velocity  field  is  investigated  that  is  possible  under  a uniform 
state  of  stress  at  the  yield  limit*  Two  cases  are  distinguished 
according  to  whether  all  principal  components  of  the  stress 
deviation  are  different  from  sero  ("regular  case")  or  not  ("degen- 
erate case")*  The  analytical  description  of  the  velocity  field 
involves  two  arbitrary  functions  in  the  degenerate  case  but  no 
such  functions  in  the  regular  case* 


1*  Introduction*  The  theory  of  the  slip  line  field 
in  plane  plastic  flow  is  one  of  the  best  developed  branches  of 
the  mathematical  theory  of  perfectly  plastic  solids*  Among  the 
problems  that  have  been  solved  in  this  field,  those  involving 
rectilinear  boundaries  are  predominant*  In  many  of  these  prob- 
lems, the  slip  line  field  contains  finite  regions  of  uniform 
stress*  The  most  general  velocity  field  possible  in  such  a re- 
gion is  obtained  by  the  superposition  of  two  arbitrary  shear 
flews  in  the  directions  of  the  maximum  shearing  stresses*  The 
analytical  description  of  this  veloolty  field  Involves  two 
arbitrary  functions  of  one  variable  each* 

Whereas  fields  of  plane  plastic  flow  have  been  studied 
systematically  (see,  for  Instance,  Ref*  1,  Chaps*  6-9,  or  Ref*  2, 
Chaps*  5-7) f three-dimensional  plastic  flow  is  practically  unex- 
plored* Finsl  (Ref*  3)  and  Thomas  (Ref*  4)  obtained  some  general 


The  results  presented  in  this  paper  were  obtained  in  the  course 
of  research  conducted  under  Contract  N7onr-35801  between  the 
Office  of  Naval  Research  and  Brown  University* 
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results,  but  did  not  construct  examples  of  three-dimensional 
flow 8*  De  Simon!  (Ref*  ?)  attempted  to  construct  a three-dimen- 
sional plastic  flow  field,  but  closer  Inspection  of  his  example 
reveals  this  to  be  a field  of  plane  flow  referred  to  coordinate 
axes  that  are  oblique  to  the  plane  of  flew*  A non-trivial  and 
genuinely  three-dimensional  field  of  plastic  flow  is  found  in  a 
paper  by  Hill  (Ref*  6)*  This  flow  field  describes  the  incipient 
plastic  flow  in  a prismatic  bar  made  of  a plastic  rigid  material 
and  subjected  to  combined  tension,  torsion,  and  bending*  The- 
equatlona  for  the  velocity  components  (Ref*  6,  Bqs*  17)  contain 
an  unknown  function,  however,  and  the  non-linear  differential 
equation  for  this  function  has  not  yet  been  solved*  When  this 
function  is  set  equal  to  sero  and  the  necessary  adjustment  of 
constants  is  made,  Hill's  flow  field  reduces  to  a special  case 
of  the  "regular"  field  discussed  in  the  present  paper* 

In  view  of  the  role  that  regions  of  uniform  stress  play 
in  the  theory  of  plane  plastic  flow,  it  seems  worthwhile  to  inves- 
tigate the  most  general  three-dimensional  flow  field  that  is 
possible  under  a uniform  state  of  stress  at  the  yield  limit*  It 
will  be  seen  that,  oontrary  to  what  is  the  case  for  plane  flew, 
the  analytical  description  of  the  moat  general  three-dimensional 
flow  field  possible  under  uniform  stress  does  not,  in  general, 
involve  any  arbitrary  functions* 

2*  Fundamental  equations.  Choose  the  rectangular  axes 
X,  Y,  and  Z,  in  the  principal  directions  of  the  uniform  stress 
field  and  denote  the  given  principal  components  of  the  stress 
deviation  by  K,  L,  and  M,  and  the  unknown  velocity  oomponents  by 
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U,  V,  and  W*  According  to  Mlses*  theory  of  plasticity  (Ref*  7)* 
the  components  of  the  velocity  strain  at  a generic  point  are 
proportional  to  the  components  of  the  stress  deviation  at  this 
pointy  the  non-negative  factor  of  proportionality  being  a func- 
tion of  position*  Applied  to  the  present  oaae9  this  flew  rule 
leads  to  the  following  equations  in  vhioh  subscripts  denote  dif- 
ferentiation with  respect  to  the  coordinates! 

*i  XK,  Vy  «=  XL,  Wz  * XHj  (1) 

Uy  + Vx  = 0,  V2  4 WY  * 0,  Wx  + Uz  = 0*.  (2) 

In  these  equations,  X denotes  the  non-negative  factor  of  propor- 
tionality* 

Since  K,  L,  and  M,  are  constants,  elimination  of  the 
velocity  components  between  Bis*  (1)  and  (2)  leads  to  the  follow- 
ing differential  equations  for  the  unknown  function  Xs 

KX^y  = 0,  1X22  * 0,  MXgg  + KX22  c 0}(3) 

KXy2  * 0,  IiXyy  ® 0,  MXjy  B 0#  (4) 

These  are  the  equations  of  compatibility  (see,  for  instance.  Ref* 
8,  p*  27,  Bis*  10*10)  for  the  components  of  the  velocity  strain 
when  these  are  expressed  in  terms  of  the  components  of  the  stress 
deviation  and  the  factor  X* 

In  discussing  Eqs*  (3)  and  (4),  we  must  keep  in  mind 
that  K,  L,  and  H,  are  the  principal  components  of  the  stress 
deviation*  Accordingly, 


K + L + M « 0* 


(5) 


The  following  two  cases  must  he  discussed  separately. 

1)  None  of  the  quantities  K,  L,  and  M,  vanishes*  Equations 
(3)  then  yield 


B ^YY  r'  ^22  = °* 

(6) 

2) 

This  will  be  called  the  regular  case* 

One  of  the  quantities  K,  L,  and  M,  vanishes* 
instance,  H = 0,  we  have  L = - K,  by  (5),  and 
and  (4)  furnish 

If,  for 
Eqs.  (3) 

XXX  " *YY  ~ °*  X22  = 0; 

(7) 

«r 

N 

II 

O 

— 

>* 

s 

II 

o 
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(8) 

This  will  be  called  the  degenerate  case^ 

For  the  regular  case,  Eqs*  (4)  and  (5)  yield 

XcAX  + BY+CZ  + D, 

(9) 

where  A, 
case,  it 

B,  C,  and  D,  are  arbitrary  constants*  For  the  singular 
follows  from  Bqs*  C7)  and  (8)  that  X has  the  form 

X » f(X  ♦ Y)  ♦ g(X  - Y)  + C 2, 

CIO) 

where  C again  denotes  an  arbitrary  constant* 

3*  Regular  velocity  fields.  With  X as  given  by  Eq* 

(9),  Integration  of  Eqs*  Cl)  and  (2)  will  yield  a velocity  field 
that  Is  compatible  with  the  considered  uniform  stress  field*  Since 
X is  non-negative,  hewever,  such  a velocity  field  can  be  construct* 
ed  only  In  the  half*space  in  which  the  right-hand  side  of  (9)  Is 


non-negative*  The  following  analysis  is  simplified  by  the  intro- 
duction of  new  rectangular  coordinates  x,  y,  and  st  that  are 
chosen  so  that  this  half- space  corresponds  to  a 2 0* 

Let  ky  -t,  and  m,  and  p,  q,  and  r,  be  the  normal  and 
shear  components  of  the  stress  deviation  with  respect  to  the  new 
coordinate  axes,  and  u,  v,  and  v,  the  components  of  the  velocity* 
Equations  (1)  and  (2)  must  now  be  replaced  by 


«x 

= Xk, 

▼y  ' 

= *»» 

(11) 

Uy  + Vx 

* 2Xp, 

+ Wy  = 2X1, 

w_  + u = 2Xr, 

A E 

(12) 

where 

X = CE,  (c  > 0). 

(13) 

To  within  a velocity  field  that  corresponds  to  a rigid 
body  motion,  the  most  general  velocity  field  compatible  with 
thes*  equations  is  given  by 


u = o(kxz  + pys  + rs;, 


v = 


p 2 

c(pxz  + ly z + qz  ) 

- £ c(kx2  + -ty2  + 2pxy 


r (*  > 0) 


2 

ms  )• 


(l»f) 


This  velocity  field  involves  no  arbitrary  functions,  in  fact  no 
arbitrary  constants  other  than  o* 

Por  given  values  of  the  components  of  the  stress  devia- 
tion, the  velocity  field  (14)  is  valid  only  in  z 2 Under  cer- 
tain conditions,  it  is  possible,  however,  to  continue  this  field 
into  z < 0 by  admitting  different  uniform  states  of  stress  on  the 
two  sides  of  the  plane  z = 0,  We  use  the  prime  to  distinguish 
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quantities  in  t < 0 from  the  corresponding  quantities  in  z > 0. 
Since  X cannot  be  negative,  the  constant  c*  must  be  negative*  If 
the  velocity  field  (14)  is  to  be  valid  throughout  space,  we  must 
therefore  set  k1  = - k,  V = - -t,  •••  , r*  = -r*  Equilibrium 
between  the  tractions  transmitted  across  s = 0 requires,  however, 
that  q*  = q,  rf  = r.  There  is  no  such  condition  for  m and  m* 
because  the  normal  stresses  in  the  z direction  are  obtained  by 
adding,  to  m and  m*,  the  respective  mean  normal  stresses  which 
do  not  appear  elsewhere  in  our  analysis*  It  follows  from  the 
preceding  discussion  that 

q s q * = r = r*  =0*  (15) 

With  this  restriction,  the  velocity  field  (14)  can  be  continued 
into  z < 0*  The  plane  z = 0 is  then  a discontinuity  surface  of 
the  stress  field  (for  a discussion  of  such  discontinuity  surfaces 
see  Ref*  9)«  The  application  of  this  type  of  stress  and  velocity 
fields  to  a plastic  plate  under  uniform  bending  and  twisting 
moments  is  obvious* 

4*  j&ggaasfl&fi. Jg&lflfil ta-X With  \ as  given  by 
Bq*  (10),  integration  of  Eqs*  (1)  and  (2)  yields  the  following 
result:  to  within  a velocity  field  that  corresponds  to  a rigid 
body  motion,  the  most  general  velocity  field  possible  in  the  de- 
generate case  is  given  by 

U = K [F(X  + Y)  + G(X  - Y)  + CXZ], 

V = -K[F(X  + Y)  - G(X  - Y)  + CYZ],  (16) 

W = - ^ CK(X?  - Y2), 
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where  the  functions  F and  0 have  as  derivatives  the  functions  f 
and  g appearing  in  (10)#  This  velocity  field  is,  of  course, 
restricted  to  the  region  in  which  \ as  given  by  Eq.  (10)  is  non- 
negative* 

The  known  results  for  plane  plastic  flow  under  constant 
stress  are  obtained  from  (16)  by  setting  C = 0*  The  general 
degenerate  field  (16)  results  from  the  superposition  of  the  field 
of  plane  plastic  flow  represented  by  the  functions  F and  G,  and 
the  field  represented  by  the  terms  with  C in  Eqs«  (16).  The 
velocities  of  the  latter  field  tend  to  tilt  the  planes  X = const, 
or  Y = const,  about  their  intersections  with  the  plane  Z = 0 and 
to  transform  the  planes  Z = const,  into  congruent  hyperbolic 
paraboloids. 
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New  York  University 

Institute  for  Mathematics  & Mechanics 

45  Fourth  Avenue 

New  York  3*  New  York 

Attn:  Professor  R.  C our ant 

Dr*  0.  Hudson 

New  York  University 
New  York  53  > New  York 

Attn:  Dr.  C.  T.  Wang,  Department  of  Aeronautics 

Northwestern  University 
Evanston,  Illinois 

Attn:  Dr*  M.  M.  Hetenyl 

University  of  Notre  Dame 
Notre  Dame,  Indiana 
Attn:  Dr*  P.  A.  Beck 

Ohio  State  University 
Columbus,  Ohio 

Attn:  Dr.  B.  A.  Boley 

Pennsylvania  State  College 
State  College.  Pennsylvania 
Attn:  Dr*  M.  Gensamer 

Dr.  J*  A.  Sauer 
Dr.  Joseph  Marin 
Dr.  J.  W.  Fredrickson 

Princeton  University 
Princeton,  New  Jersey 
Attn:  Dr*  S.  Lefschets 

Dr*  L.  Lees 
Dr*  J.  V.  Chary k 

Rensselaer  Polytechnic  Institute 
Troy,  New  York 
Attn:  Library 

Dr*  Paul  Leiber 

Santa  Clara  University 
Santa  Clara,  California 
Attn:  Dr*  R*  M*  Hermes 

Stanford  University 
Stanford,  California 
Attn:  Dr*  L*  Jacobsen 

Dr*  A.  Phillips,  Department  of  Mechanical  Engineering 
Dr*  J.  N.  Goodler 
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Stavena  Institute  of  Technology 
H&boken,  New  Jersey 
Attn*  Dr*  E.  0.  Schneider 

Swarthmore  College 
Swarthmore,  Pennsylvania 
Attn*  Capt.  W.  P.  Roop 

University  of  Texas 
Austin  12,  Texas 

Attn*  Dr.  A.  A.  Topractsoglou 

University  of  Utah 
Salt  Lake  City,  Utah 
Attn*  Dr.  H.  Eyring 

Washington  State  College 
Pullman,  Washington 
Attn*  Dr.  B.  Pried 

Wheaton  College 
Norton,  rir.ssachusetts 
Attn*  Dr.  H.  Geiringer 

Aerojet,  Inc* 

Asusa,  California 
Attn*  7.  Zwicky 

Aluminum  Company  of  America 
New  Kensington,  Pennsylvania 
Attn*  R.  L.  Templln 

Armstrong  Cork  Company 
Lancaster,  Pennsylvania 
Attn*  J.  W.  Scott 

Bell  Telephone  Laboratories 
Murray  Hill,  New  Jersey 
Attn*  C.  Herring 

J.  M.  Richardson 
D.  P.  Ling 
W.  P.  Mason 

Corning  Glass  Company 
Corning,  New  York 
Attn*  J.  T.  Littleton 

E*  I.  Dupont  de  Nemours  & Co* , Inc. 

Wilmington  98,  Delaware 

Attn*  J.  H.  Paupel,  Materials  of  Construction  Section 

General  Electric  Company 
Schenectady,  New  York 
Attn*  H.  Fehr 

H.  Poritsky 
7.  H.  Hollomon 


7 


General  Motors 
Detroit,  Michigan 
R-l  Attns  J.  0.  Almen 

Lockheed  Aircraft  Company 
Department  ?2-25,  Factory  A-lt  Building  66 
Burbank,  California 
R-l  Attn:  Engineering  Library 

Midwest  Research  Institute 
Kansan  City,  Missouri 
R-l  Attn:  C.  0*  Dohrenwerxd 

R-l  M.  Oolan 

Pratt  d Whitney  Aircraft  Corporation 
East  Hartford,  Connecticut 
R-l  Attn:  R*  Morrison 

U.  S*  Rubber  Company 
Passaic,  New  Jersey 
R-l  Attn:  H.  Smallwood 

Welding  Research  Council 
Engineering  Foundation 
29  West  39th  Street 
New  York  18,  New  York 
M-l  Attn:  W.  Spraragen,  Director 

Westinghouse  Research  Laboratories 
East  Pittsburgh,  Pennsylvania 
M-l  Attn:  Dr.  A.  Nadai 

R-l  Dr.  E.  A.  Davis 

Westinghouse  Electric  Corporation 
Lester  Branch  P.  0. 

Philadelphia,  Pennsylvania 

R-l  Attn:  R.  P.  Kroon,  Mgr.  of  Engineering,  AGT  Division 
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